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lO ■ Abstract 

o . 

^\ . New classical integrable systems of Camassa-Holm peakon type are proposed. They realize 

the maximal even piecewise--D2 generalization of the Calogero-Frangoise flows, yielding pe- 
riodic and pseudoperiodic trigonometric/hyperbolic potentials. The associated r-matrices 
^ \ are computed. They are dynamical and depend on both sets {pi} and {qi} of canonical 

' variables. 
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1 Introduction 

Classical r-matrices for the Lax representation of Liouville-integrable systems, depending on 
some of the canonical variables (so called "dynamical r-matrices"), have been the subject of 
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a very extensive set of studies in these past years, starting with the first mention of their 
relevance [1], and their use in non- ultra local Poisson algebras [2, 3]. 

Most studies afterwards have been concerned with a subset of such matrices, obeying the 
so-called "dynamical Yang Baxter algebra" . A paradigmatic example occurred in the Calogero- 
Moser r- matrix formulation [4] , and classification theorems were extensively developed, see for 
instance [5] . These matrices depend only on one half of Poisson- commuting canonical variables 
(either momenta or positions), identified in the general theory [6] as coordinates on the dual of 
an Abelian Cartan algebra. The non-abelian situation is discussed in e.g. [7]. 

The above-mentioned associated dynamical classical Yang Baxter equation, obtained as a 
sufficient consistency condition for associativity of the Poisson brackets, can then be related 
to dynamical quantum algebras such as the Gervais Neveu Felder (GNF) algebra [8, 6] or the 
dynamical reflection algebras [9, 10, 11]. 

More precisely: the clearest case is provided by the Ruijsenaars- Schneider quadratic Poisson 
structure involving two matrices [12, 13]. One is then able [13, 10, 9] to identify the consis- 
tency equations as a classical limit of the Yang-Baxter equations for a dynamical reflection 
algebra. They include in particular the quantum Gervais Neveu Felder cubic equation (sym- 
bolically written RR^R = R^RR^) for e.g. Interaction Round a Face (IRF)-type quantum 
i?-matrices, and its associated adjoint equation RSS^ = SS^R. It is interesting to remark that 
the flrst mentioned example, i.e. the scalar Calogero Moser r-matrix, obeys in fact a consistent 
combination of these two equations. 

Even in such cases when the r-matrix or matrices depend on Poisson-commuting dynamical 
variables, the precise derivation of the associated classical Yang Baxter equation is thus a 
subtle procedure requiring the explicit knowledge of the Lax matrix and r-matrix. Obtaining 
the specific form associated with the GNF equation is not a priori guaranteed. 

By contrast there are fewer results regarding r-matriccs depending on non-Poisson commut- 
ing sets of dynamical variables (hereafter loosely denoted "p, q dependent r-matrices" ) . Indeed 
few examples are known, even less explicitely so, no associated classical "dynamical" Yang 
Baxter equation has been derived, let alone any quantization thereof. Note at this point that 
we expect any such quantization to have very unusual features compared to the "standard" 
dynamical case, involving as it should non- commutative deformation parameters p, q (instead 
of the sole c-number parameters q) inside the generalized R matrix. 

The study of such r-matrices is thus clearly a very challenging and promising issue in 
the domain of integrable systems and this paper will present an attempt at expanding the 
current state of knowledge about them by constructing new examples of Lax matrices and their 
associated r-matrices. 

To position the problem, let us flrst recall a few known examples where the r-matrix exhibits 
these dependance: the complex sine Gordon case [3] (explicit pair r, s but no decoupled Yang 
Baxter equation); the elliptic Galogero Moser model [14, 15] in its spectral parameter- free 
formulation [16] (nonexplicit); the BCn Ruijsenaars Schneider Lax matrix [17] (nonexplicit) ; 
and possibly the recently proposed 2x2 Lax matrix associated with special solutions of the 
Gamassa Holm equation [18]. 

The simplest example (albeit nonexplicit yet) is provided by the Ragnisco-Bruschi [19] 
construction of the r-matrix for the Lax matrix yielding the Calogero-FranQoise (CF) Hamilto- 
nian [20]. The Lax matrix reads: 

L = Ljj-eij ; = y/p^ A{qi - qj) ; with A{q) = cosh(i/q/2) psinh(i/ |q/2|) (1.1) 
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The CF Hamiltonian is Tr(L^). This dynamical system may be viewed as a generahzation 
of both Camassa-Holm peakon integrable system (where A{z) — exp{iy\z\)), describing the 
dynamics of specific singular "peaked" solutions to the famous Camassa Holm shallow wave 
fiuid equation [21, 22, 23], and the related Calogero system [24] (where A{z) — cos{z)). 

We shall be classifying here the Lax matrices a la Calogero- Prangoise admitting an r-matrix. 
We shall in this paper restrict ourselves to Lax matrices where the diagonal potential An is a 
constant (we shall comment upon this at the end) and A is an even function. 

The problem will then be solved by first imposing and solving a necessary condition on 
A derived from the vanishing of the Poisson bracket {TrL,TrL^} and {TrL^,TrL^}; then 
explicitely constructing the r-matrix will show that any such solution be indeed integrable. 
Some further comments and perspectives shall be added at the end. Not surprisingly the 
periodic peakon potential [25] will be obtained as one solution. 



2 The model 

We assume the following form for the Lax matrix: 



E, . . ^ \ Aii — Qi some complex constant ,„ ^ . 

VMJA,e;, w,th^ ^,;, = ^(,._,,)fori^j (2.1) 



As usual eij denotes the n x n matrix representation of gl{n). We require that A be an even 
complex function of a real variable, with isolated singularities (more precisely separated from 
any other such point by a distance strictly bounded from below), at least Di around regular 
points and such that be a piecewise D2 function. 



2.1 Necessary conditions 

Starting from the formulae for the first three Hamiltonians defined from L: 

i i 

H2 = TrL^ = Y^PiPi Hi = Y.pI + Y^PiPj ^\<li - Qj) (2-2) 
H3 = TrL^ = ^PiPjPk Aij Aj^k Ak,i = 
^Pi <4 + 3 ^pfpj ai A^{qi - qj) + ^ PiPjPk A{q^ - q^ A{qj - qk) A{qk - qi) 

we impose simultaneous vanishing of all Poisson brackets between them as a necessary 
condition for integrability. 

With the convention {pj, qj} — 5ij, one gets 

{ifi, H2} = 2 ^PiPj {ai - ttj) A{qi - qj) A'{qi - qj) 

Hi and H2 are in involution if and only if either the potential is a piecewise constant function (a 
trivial solution, so we will solve excluding this case) or the diagonal part of the potential is uni- 
form: Oj = a. Setting all Oj's to a the Poisson bracket of H2 and H^, after some rearrangements. 



3 



reduces to 

{H2, Hs}^6j2 PiPj PkPi ^k,i {Ai,i Aj,k - Ajj {Ai,i A[^ - Aij A[,) 
with A'- J = A'{qi — qj) for i ^ j and A'- - = . In particular the p^PhPc coefficient is a multiple 

of: {A,,l Aa,c - a A5,e) [{Alb - <c) 4,c + K,b + Aa,c K^,) A,,, - a A^, + Aa,c ■ 

Let us denote the differences in positions in terms oi x — Qa — Qb and y — Qc — Qa', the 
commutation of H2 and yields: 

A(x)A(y)^aA(y + x) (2.3) 
or 

(A^ix) -A\y)) A'iy+x) + (AA'iy) -AA'{x)) Aiy+x) = a (^(x) A'iy) -Aiy) A'ix)) (2.4) 

Note that, if for x and y in the vicinity of some Xo and yo equation (2.3) is satisfied, then 
A'{x)A{y) = A{x)A'{y) = aA'{y + x) yielding also A"^ (x) A' {y + x) = AA'{x)A{y + x) and 
A^{y)A'{y + x) = AA' {y) A{y + x) . So equation (2.4) has to be satisfied in the full plane x,y 
except perhaps in some isolated points (necessarily singular for A). 

Introducing the variables s — y + x and d — y — x, as well as the three functions 
D{s, d) = A^{y)-A^{x), F{s, d) ^2A{y) A{x) and G = daF, equation (2.4) now reads: 

dsD{s, d) A{s) - D{s, d) A'{s) = a G{s, d) (2.5) 

In the vicinity of a regular point, the derivation with respect to d also yields: 

dddsD{s, d) A{s) - ddD{s, d) A\s) = a ddG{s, d) (2.6) 

Eliminating A'{s) from equations (2.5) and (2.6) and defining A = DdddgD — ddDdgD, one 
gets: 

A(s, d) A{s) = a {D{s, d) ddG{s, d) - ddD{s, d) G{s, d)) (2.7) 

Note that D{s,d) = D{d,s) directly from its definition since the exchange of s and d 
corresponds to the transformation of x into —a; and A is an even function. Hence A is also a 
symmetric function and one gets from (2.7): 

A(s, d) A{d) = a {D{s, d) dsH{s, d) - dsD{s, d) H{s, d)) (2.8) 

with H — dgF (note that F is also symmetric in s and d). 

Actually, by direct computation, one obtains that the two rhs of equations (2.7) and (2.8) 
are equal. Let us denote them by aK with K the following symmetric function: 

K{s, d) = D ddddF - ddD ddF = D d^d^F - d^D d^F (2.9) 

Since we have excluded from the begining any piecewise constant potential, we obtain A = 
as a necessary condition: 

A = D dddsD - ddD d^D = (2.10) 
This very equation can be solved on any regular domain of the plane, so we get piecewisely: 

D{s,d) = a{s)(3{d) (2.11) 
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with a and /? two Di-functions on the interval of s, respectively d, where they build D. Note 
that if we change the domain of the plane where we solve the equation the functions a and /3 
are not necessarily the same functions as in the first domain, even if they happen to be defined 
at some common value of s {or d). 

Remark that since D{s, d) — A^{y) —A^{x) and since has been supposed piecewisely D2, 
a and (3 are actually D2, one thus gets piecewisely: 

\ {{Ay{y) - {A^ix)) = a"{s) /3(d) = a{s) ^"{d) (2.12) 

So for any pair of functions a and /3, there exists some complex constant v such that 

a"{s) - v^a{s) , (3"{d) = f3{d) and {A")" (y) - 4 A^y) = {AY (x) - 4. A" (x) (2.13) 

Note that the last equality holds for any couple of regular points. One finally obtains the 
following form for the square of the potential: 

A2(x) =7 + A^e2''-^ + /x^e-2'^-^ (2.14) 

with 7 some complex constant and v^, and /i^ three piecewise constant functions (to dis- 
criminate between A and /i we will assume Re{vx) > 0). Reinserting this form into (2.10) 
yields: 

{ul - ul) (A. - - + ^) (A, ^ + e'^ ^) + 4 {ul A, - ul A. /..) = (2.15) 

One obtains that is actually a complex constant, i.e. — Vy, that we will accordingly 
denote by u in the following , we can assume that it is not zero since we have already discarded 
the piecewise constant potential. More interestingly the product A-^ /i^ is also constant and we 
shall now call it (up to a normalization factor). 

One then easily shows that p = is essentially not compatible with the assumption of an 
even potential and leads to A constant. In the following we will thus assume p 7^ 0. 

2 

Let us now introduce another piecewise constant function 9^, parametrize Xx fJ^x = fg by 
Ac. = f e-"^- and /x^ = f e'^''-, and get: 

A'^{x) = p (cosh2(i/ {x - 9x/2)) - cosh^{v Zo)) 
= p sinh(i/ {x + zq- Bx/2)) sinh(i/ {x - zq - 9x/2)) (2.16) 

with Zo a complex constant. 
This expression of yields: 

D{s, d)^p (cosh2(i/ (y - ey/2)) - cos\?{v {x - Ox/^))) 
^ psmh(u{s- (ey + ex)/2))sinh(u(d- (ey-ex)/2)) (2.17) 

Using this expression in equation (2.5), one shows that: 

dd{A{x) A{y)) = A{s) usmh{u {d - {9y - ^,)/2)) 

with A{s) some function of s. Hence there exists some function i3 of s such that: 

A(x) A(y) = A(s) cosh(i/ (d - {Oy - e^) /2)) + B{s) (2.18) 
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From (2.16), one also gets: 

A'ix) A'iy) = p2 sinh(^ {x + zo- 9J2)) smh{u {y + Zo - 9y/2)) 
smh(i/ {x — zq — 9x/2)) smh(i/ {y — zq — Oy/2)) 
= pV4 (cosh(i/ (s + 2 z„ - {Oy + ^,)/2)) - cosh(i/ {d - {Oy - ^,)/2))) 

(cosh(i/ (s - 2 - (^j, + e,)/2)) - cosh(i/ (d - (^^ - e,)/2))) (2.19) 

The compatibihty between (2.18) and (2.19) imposes: 

cosh(i/ {s + 2zo- {9y + 9x) /2)) = cosh(i/ (s - 2 - (^j, + e^) /2)) 

or equivalently sinh(z/2zo) = restricting the constant cosh^(z/Zo) to be either or 1. Up to a 
shift on ^j., we have obtained that the potential is of the following form: 

A{x) = Ao cosh(i/ {x - 9x12)) (2.20) 

with Ao a complex constant. Equation (2.5) now reduces to an equation on the piecewise 
constant function 0^ (defined modulo AiTi/v): 

Ao cos\i{v {9 x + 9y-ex+y)/ 2))^ a (2.21) 

The only other condition is the parity of A, imposing d-x = —9x- In particular this is true for 
a; = where A is thus continuous, but generically not differentiable. 

So either 9x — = —9-x for all positive x, or it admits a first strictly positive discon- 
tinuity in xi, and reads 9x — (2^2,-1-1) ^o+, with Ux the integer quotient of x divided by xi 
{x ^UxXi^ Tx, Tx e [0, xi[ and Ux G Z). 



Finally, we have established the following proposition: 



Proposition 1 

A Lax matrix of type (2.1) which defines a completly integrable system can at most depend on 
three complex parameters {A^,!/ and ^o+) and possibly a strictly positive real one {xi) It has 
to be of the following form: 



with either 9x — sign{x) or 



XI VM" AoCosh{v {qi -qj-\ 9^q^-q.))) ey, (2.22) 

Ox = (2 + 1) 6*0+ for x = n^xi + rx, 
Tx G [0,a;i[, and Ux G Z 



The first set of solutions with 9x — sign{x) 6*0+ corresponds exactly to the Calogero-Frangoise 
Lax matrix (1.1). Note that within the new set of solutions, if one put xi — ^o+, one gets a 
non trivial continuous periodic potential. 

2.2 The r-matrix 

We will now estabhsh Liouville-integrability for the Hamiltonians associated with the Lax 
matrix (2.22) by explicitely constructing the r-matrix. 
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Following closely the r-matrix derivation in [19] it is easily seen that it hangs upon the 
following technical lemma: 

Lemma 

If there exist two odd functions B and C, such that the even function A obey the consistency 
equation (labeled (2.10) in [19]) where again s — x + y : 

ds{A{x) A{y)) = Bis) + Ais) {C{x) + C{y)) (2.23) 

the Lax matrix L = \/PiPj ^{Qi ~ Qj) ^ij has a linear r-matrix Poisson structure: 

{Li,L2} = [ri2, Li] — [r2i, L2] with the standard notations Li = L (g) 1 and L2 = 1 <8) L 

and ri2 = ^ rijki Cy O Cki ; = ^ rijki O given by: 

ri2 = ^ C{qi - Qj) Bij (eij + eji) + ^ '^^'^ ^^'J ~ ^J'^ ® ^^'^ + ^^-^"^^ 

Sjf^ are the elements of any symmetric matrix S solving the algebraic equation: 

[S, L]ij = .Jp'p'j B{qi - Qj) (2.25) 



Proof: Proof can be obtained from direct verification, once it is established that the Pois- 
son brackets of L with itself are given uniquely in terms of the function ds{A[x) A{y)) . 

Remcirk 1 

Existence of a solution S to the algebraic equation (2.25) is obtained by rewriting the defin- 
ing equation in an orthogonal basis for the symmetric matrix L when L is real. If L is complex 
symmetric, one needs to assume separately that it can be diagonalized (which is true if L is in 
general position). In both cases one also needs to assume that the eigenvalues are non degen- 
erate. Note that the issue of degenerate eigenvalues for a Lax matrix is in any case a delicate 
one and needs a case-by-case treatment. 

Remcirk 2 

S is in any case not unique, since in particular any polynomial in L yields a symmetric 
matrix commuting with L itself. 

We are currently working on an explicit formulation of solutions to equation (2.25). 

We now state the key result of this section: 

Proposition 2 

The function A{x) in (2.20) obeys the r-matrix consistency equation (2.23) with the functions: 

B{x) = i(A'(x) A{0) - ^ A{x) A'iO)) (2.26) 
= ^(sinh(i/ (x - e^/2)) cosh(i/^o+/2) + ^ cosh(i/ {x - 9^/2)) smh{v ^0+72)) 
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and ^(^) = ^^^'(0) = -^^sinh(^^o,/2) (2.27) 

Oddness of functions B and C follows from = —Ox- The proof proceeds by direct 
checking of the "addition" formula (2.23), paying due attention to the consistency of this 
addition for the step function O^- 

Hence all necessary potentials A yield a Lax matrix which exhibits an r-matrix structure of 
the form obtained by Ragnisco-Bruschi, albeit with periodic or more generally pseudoperiodic 
coefficient functions B and C. We have thus completely classified the integrable classical Lax 
operators of the Calogero-Frangoise form. 

3 Remarks and perspectives 

1. We have obtained a new set of pseudoperiodic potentials as a result of our systematic 
search for even Calogero-FrauQoise type potentials. It would be interesting to know how to 
formulate alternatively algebraic integrability properties (extension of the Abel maps and asso- 
ciated hyperelliptic curves [25]) for such pseudoperiodic potentials. In principle the conserved 
Hamiltonians should be moduli of the associated curve. The possible connection with the 
pseudoperiodic solutions derived in e.g. [23] may also be interesting. 

Of course among these potentials one recovers the fully periodic peakon potential derived 
by Beals et al [25]. Here we prove its r-matrix (hence Liouville) integrability whereas [25] had 
proven algebraic integrability in the above sense. 

2. Interpretation of the pseudoperiodic "pcakons" in the frame of Camassa Holm equation 
may be interesting, possibly as solutions for a system in a finite volume with some specific 
boundary conditions. 

3. It is easily checked that a folding procedure exists which consistently yields new integrable 
systems from these initial pseudoperiodic ones. Starting from the N — 2n(-|-l)-body system 
with positions labeled now from —n to -\-n (without/with 0) one proves immediately that 
the folding constraints pi — p-i, qi = q^^i are invariant by evolution under all even (resp. 
odd) degree Hamiltonians. Hence the reduced dynamical system itself has n invariant Poisson 
commuting Hamiltonians under the reduced Poisson bracket. Its Lax matrix takes the form: 

n 

L = ^ Lij (ejj + e_ij + _j + e_i _j) with: 

i,j=m 

L,, = V^(cosh^^^^y^ + psinh^^^i^) (3.1) 

Here one first (rather trivial) possibihty of extension of the strict Calogero Prangoise frame- 
work is made explicit. 

4. In order to proceed towards a deduction of the form for the Yang-Baxter type equation 
one key point now is to get an explicit form for the r-matrix by solving the matrix equation for 
S. Partial results are now available and we hope to soon have the general form of this r-matrix. 
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